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VO ' Divergent carrier-density fluctuations equivalent to the critical opalescence of gas-liquid 

transition emerge around a metal-insulator critical point at a finite temperature. In con- 
trast to the gas-liquid transitions, however, the critical temperatures can be lowered to zero, 
which offers a challenging quantum phase transition. We present a microscopic description 
of such quantum critical phenomena in two dimensions. The conventional scheme of phase 
transitions by Ginzburg, Landau and Wilson is violated because of its topological nature. 
It offers a clear insight into the criticalities of metal-insulator transitions (MIT) associated 
, - with Mott or charge-order transitions. Fermi degeneracy involving the diverging density fluc- 

r^ . tuations generates emergent phenomena near the endpoint of the first-order MIT and must 






o 



shed new light on remarkable phenomena found in correlated metals such as unconventional 

cuprate superconductors. It indeed accounts for the otherwise puzzling criticality of the Mott 

CN . transition recently discovered in an organic conductor. We propose to accurately measure 

>■ . 
I — . enhanced dielectric fluctuations at small wave numbers. 
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Insulating phases appear for various reasons, ^"^ ranging from simple band insulators, An- 
derson insulators caused by randomness, to insulators driven by Coulomb repulsion of electrons 
'^ ' such as Mott insulators. The band and Anderson insulators always emerge through continu- 

Q I ous MITs from metals at temperature T = 0, meaning the zero critical temperature Tc = 0. 

However, Mott transitions and charge-order transitions often show first-order character with 

huge jumps in conductivity.^ The jump decreases with increasing temperature and terminates 
/S ■ 

;_( ■ at the critical point at T = Tc > 0. In this letter, we explore the mechanism for the emergence 

C^ ' 

of such first-order and continuous MITs. 

In the simplest scheme by Ginzburg, Landau and Wilson (GLW) for phase transitions,'^'^ 

the free energy F on the mean-field level is expressed as a function of the order parameter m 

as 

F'i^) = ^^GLm^ + jBcLm^ - mh, (1) 

where we omit the spatial dependence. Here, Bqi is regarded as a positive constant, whereas 
a phase transition occurs through the sign change of Aql, which depends on the control 
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parameter g around the critical point gc as Aql = AcLoig — gc)- Within the mean- field theory, 
when the external field changes the sign from h > to h < 0, the free-energy minimum point 
jumps at /i = from m-^. = niQ to m_ = —mo for g < gc, which determines a critical 
exponent /? defined by Am = \m^ — m_| (x \g — gd^ as (3 = 1/2. When g is kept at gc and 
h is controlled, another exponent 6 defined by m oc h^' becomes 6 = 3. The susceptibility 
X = linih^Qm/h (x \g — ^d ''' determines the third critical exponent 7 = 1. These mean-field 
values are modified under fluctuation effects while the mean-field exponents are correct above 
the upper critical dimension dc = 4. This Ising universality class also applies to the classical 
gas-liquid transition. The order parameter of the gas-liquid transitions is regarded as the 
particle density measured from the critical point. 

A recent study of (Vi_a;Cr2:)203 by Limelette et al.^ has suggested that the metal-to- 
Mott-insulator transition also belongs to the Ising universality.'^ Similarly to the gas-liquid 
transition, the order parameter is interpreted as the carrier density X, which drives the sin- 
gularity of conductivity. In analogy with the critical opalescence in the gas-liquid transition, 
diverging carrier density fluctuations control physics around the critical point. 

However, in earlier studies, a scaling theory of a different universality class has been 
proposed when Tc can be suppressed to zero,^'^ in which the simple GLW scheme is violated 
and unusual exponents 5 = 2, 7 = 1 and j3 = 1 are derived. "'^'^ A recent experiment ^-"^ on 
an organic conductor, K-(ET)2Cu(Cn)2Cl, has supported the existence of this unconventional 
universality. 

In this letter, we show that such a new scheme of quantum critical phenomena can indeed 
be constructed from a microscopic model. When the Ising-type critical fluctuations become 
involved in the Fermi degeneracy region, it deeply modifies the classical nature and a dif- 
ferent type of diverging density fiuctuations shows up with emerging effects, indicating the 
breakdown of the GLW scheme. An unforeseen region of marginally quantum critical (MQC) 
behavior sandwiched by the regions of the Tc = critical line and the first-order transitions 
may signal various instabilities and competing orders. 

To understand the essence, the Hubbard model on a Ng = L x L square lattice defined by 

ka i (ij) 

is employed with the hopping terms for the nearest neighbor ^1 = —2t{coskx + cosky) and 
the next-nearest neighbor ,^2 = 4t' cos kx cos^ . We assume < t' < t. The standard notation 
is employed their Fourier transforms are denoted by cl.^{cka)- The number operator is defined 
by fii = ^^jTiia- and Ui^ = cJ^Oia and the chemical potential is /U. The onsite and nearest- 
neighbour Coulomb repulsions are denoted by U and V, respectively. Hereafter, we take the 
energy unit as t = 1. 

We first consider the electron density at half filling, n = 1. When 81^ > U, this model 
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has a tendency of charge ordering, where electrons become alternatingly empty and doubly 
occupied at every other site. For U > 8V, one electron tends to occupy each site leading 
to the Mott insulator and the exchange interaction stabilizes the antiferromagnetic order. A 
dimensionless coupling g = (8V — U)/2t for 8V > U and g = U/2t for U > 8V is introduced 
in the following discussion. 

When one assumes the antiferromagnetic or charge orders, the Hartree-Fock approxima- 
tion (HFA) in the thermodynamic limit Ng ^ oo yields the free energy at temperature T 
represented by the order parameter A as 

nA) = -^^log{l + e-^^('y^)+^ + ^,n (3) 

with the quasiparticle dispersion 



E±{k)=Uk)± VUW+^ - f^n, (4) 

where we take the upper (lower) sign for the wavenumber satisfying ^i > 0(< 0), respectively. 
Here the order parameter A is defined by A = gm with m = j^ Si(^JT ~ ?T.jj,)e*'^'"' for the 
antiferromagnetic order in the region U > 8V and m = jj-Yli{''^i ~ 1^^)^ for the charge 
order (U < 8V), with the period being Q = (7r,7r). 

HFA ignores fluctuations and thus one might have suspicions on the reliability in strongly 
correlated systems. However, as we will see, it clarifies an unexpectedly remarkable feature 
of MIT overlooked in the literature and enables us to provide a microscopic description of 
unconventional non-GLW-type criticality. We will also discuss an estimate of the fluctuations 
beyond HFA later. 

In HFA, MIT takes place in the ordered phase, where A 7^ 0. A plays a role of opening 
a gap, but the density of states at the Fermi level remains nonzero when A < Ac = 2t', 
resulting in a metal. This is because the bottom of the upper band E^ at (vr, 0) is lower than 
the top of the lower band E^ at (7r/2,7r/2). When A increases with increasing g, MIT occurs 
at A = 2t' , above which the Fermi level lies in the gap, resulting in an insulator. Note that 
the symmetry does not change at MIT. 

We first consider T = 0. When the free energy (3) is expanded in terms of A measured 
from the transition point A = Ac, namely by y = Ac — A, we obtain 

F{Y) = F{Y = 0) + AY + ^BY^ + ^CY^ + ■■■ . (5) 

If the minimum of eq. (5) is in the region y < 0, it describes an insulator, while a metal is 
given when Y > 0. 

By using the density of states D^{E) {D-^-{E)) for the lower (upper) band, 
the coefficients in Eq. (5) are given by elementary expansion of Eq. (3) as ^ ~ 
-{l/g-Wi)Ac, B ~ -{(l + a)Q+ + (l-a)Q_} - W2 + l/g and C ~ -3^3 + 
{(l + a)(3 + a)Q'+ + (l-a)(3-a)Q'_}/Ac - (1 + afR+ + (1 - afR-, where Wi = vi, 
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W2 = vi - A^ws, and W3 = Act's - A^t;5 with 



A , , 1 



We have also introduced Q, = f ^^dE (/){£, T)Ds{E), Q'^ = J^^dE<j){E,T)Ds{Ef / Df\E), 
and Rs = /f^ dE(l){E, T)dDs{E)/dE, for s = + and - with the band cutoff A. The density of 
states for t' = is denoted by Ds (E) . Because the energy derivative of the Fermi distribution 
/ defined by (l){E,T) = -df{E)/dE = e^'^ /[T {e^^^ + l)^] is reduced to the delta function 
SiE) in the limit T ^ 0, Q^ and Q', reduce to Q^ = l)s(O) and Q'^ = DsiOf /D^^\o), 
respectively, at T = 0. Then (5+ and Q_ as well as Q'_^ and Q'_ jump from the metallic to 
insulating sides, since D^{E) and D^{E) for two-dimensional systems jump from nonzero to 
zero at the Fermi level E = Q. The factor a = —dfi/dA is defined by the relation ^ = aY 
with a = — (v f^ — t' — 2t')/{\/t'^ — t' + 2t'). The carrier density for y > is easily derived 
as 

X = 2^(1 + a)y, (7) 

which holds for t' > 0. 

Now it turns out that the coefficients B and C explicitly depend on the density of states 
at the Fermi level, which is nonzero in the metal, while it jumps to vanish in the insulator. 
Therefore, B and C jump at the MIT point at T = between Bm and Cm in the metal 
{Y > 0) and Bi and Cj in the insulator (Y < 0), respectively. We realize that such jumps of 
the coefficients violate the GLW scheme of conventional expansion. When we look into more 
details, Cj < < Cm and Bm < Bi are satisfied for < t' < t. Furthermore, within the 
realistic choice of parameters, a minimum appears at Y = Yi ^ —A/Bi < if B^ ^ ACi > 
and another minimum appears at 



Y = Ym^ {-Bm + V^m " 2ACm)/Cm > (8) 

if Bm < and B'^ — 2ACm > 0. An example of such free energy is illustrated in Fig. 1, where 
the overall form is given by connecting two cubic functions generating a double-well structure 
in total. The first-order MIT occurs when the minimum changes between Yi and Ym- 

This result indicates that electron correlation effects contain a microscopic mechanism of 
generating a first-order MIT through a synergetic and mutual feedback between the growth 
of the order parameters and MIT at the Fermi level. Although such a mechanism of the first- 
order MIT has not been well recognized in the literature, ^^ the present analytical results for 
small Y are indeed confirmed in our numerical Hartree-Fock solution of the Hubbard model, 
where the first-order MIT appears for t^g/* = 0.0557 < t'/t < 0.3645 = t'^^/t (see Fig. 2). 

With elevating temperatures, the first-order jump in the region t'^2 < t' < t^i vanishes 
at Tc^ which forms the critical line illustrated as a green curve in Fig. 2. The critical line is 
given by the crossing line of the A = surface and the Bm = surface in the parameter 
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Fig. 1. (color): Free energy of extended Hubbard model as a function of F = Ac — A. The red (blue) 
cubic curve representing the insulating region Y < (metallic region y > 0) in total form a 
double- well structure as a combined solid curve. 

space of T, g and t' /t, where the left and right minima in Fig.l merge at 1" = 0. On the 
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Fig. 2. Phase diagram of metal and insulator in parameter space of temperature T, g and t' for 
Hubbard model on square lattice. The metal-insulator boundary set by the first-order MIT is 
illustrated by the pink surface surrounded by critical line shown by green curve. The endpoints 
of the critical line are marginal quantum critical points (MQCPs) depicted by yellow circles. 
Continuous and first-order MITs at T = are shown by blue and red curves, respectively. 

other hand, in the region t' > t'^^ or t' < t^g) ^m becomes positive, where the first-order MIT 
does not exist any more, while the continuous MIT boundary further continues at T = 
as blue curves, which is given from the crossing line of two surfaces, A = and T = 0. At 
r = 0, the boundaries between the quantum critical lines (blue curves) and the first-order 
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region (pink surface) appear at t'^^ and t'^21 which are illustrated as yellow circles in Fig. 2. We 
call these two points the marginally quantum critical points (MQCPs), because they appear 
at the border (or margin) from the quantum critical line to the finite-temperature classical 
critical point. In other words, each MQCP is determined from the cross point of all the three 
surfaces A = 0, Bm = and T = 0. It is remarkable that two MQCPs appear and Tc forms a 
dome structure sandwiched by these MQCPs. The overall phase diagram illustrated in Fig. 2 
proves that a simple two-dimensional Hubbard model allows us to study the diversity of MITs 
varying from (i) the first-order, (ii) marginally quantum to (iii) continuous (quantum critical) 
ones at T = 0."*^^ 

Around MQCP, we may expand as Bm = Bm {g — Qc) or Bm = Bm {T — Tc) with T^ = 0, 
where Bm and Bm are constants. Because of Eqs. (7) and (8), the jump of X at the first- 
order MIT grows as X oc Ym oc 2Bm [g — gc) /Cm, with an unusual critical exponent (5 = 1. 
When a control parameter h such as pressure couples linearly to Y in the free energy, namely, 
if —Yh is added to the free energy at MQCP given hy A = Bm = 0, the free energy minimum 
occurs at CmX jl — /i = on the metallic side, leading to X oc Yu^' with the exponent 5 = 2. 
Figure 3 shows 5 = 2 at T = in the growth of carrier number as a function of g measured 
from MQCP at 5 = g^ where g — gc plays the role of h here. When one approaches MQCP 
strictly along A = 0, the susceptibility scales as x = dY/dh oc 1/Bm (x \g — gd^"' with 7 = 1. 
These show that the carrier density fluctuations diverge at MQCP either as x '^ \9 ~ 9c\~^ 
or as X (x X^^. We recall that the critical opalescence signals diverging compressibility in 
liquid-gas transitions. Such divergence also emerges in an electronic system even at T = 0, 
but with different exponents. All the conventional scaling relations such as (3(6 — 1) = 7 and 
the hyperscaling relation are still satisfied. 

The unusual exponents 6 = 2, (3 = 1 and 7 = 1 at MQCP are completely consistent 
with the scaling theory^ and the quantum Monte Carlo results for the filling-control MIT.^ 
Moreover they are indeed in agreement with recent experimental results obtained for a quasi- 
two-dimensional K-ET-type organic conductor. ^^ 

When Tc becomes nonzero away from MQCP, we expect the emergence of a crossover 
between the marginal quantum criticality and the conventional Ising criticalities. At T > 0, the 
jumps in B and C immediately smear out leading to effective Y dependences of B and C and 
the GLW expansion (1) ultimately becomes valid. However, at sufficiently low temperatures, 
the change from Bi to Bm, as well as from Cj to Cm, occurs within a very small interval, 
|y| < T. It turns out that if |y| > T is satisfied, the marginal quantum criticality dominates, 
while the crossover to the Ising criticality occurs around |y| ~ T. This crossover at finite 
temperatures is indeed seen in our numerical solution as seen in Fig. 3. The left panel of Fig. 4 
illustrates how Tc grows as a function of gdT = 0) — g and in the right panel we show how 
the Ising region with 6 = 3 around T^ becomes extended with increasing Tc while most of the 
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region is governed by MQCP with 6 = 2. The Ising region is very limited {\g—gc{Tc)\/t < 10^^) 
even for Tc = O.Olt, whereas the critical region governed by MQCP is wide. 

In the organic coductor/^ the exponents are measured around a finite-temperature critical 
point. However, the present unconventional criticality at MQCP may indeed be observed, 
because the experimental resolution is 0.1 MPa, while for a realistic choice of parameters 
and Tc ~ O.Olt, the Ising region appears only for the pressure P measured from the critical 
value |P — Pc\ < 10^'^ MPa, as derived from 1 MPa corresponding to change in 0.002t in 
the case of k-(ET)2[N(CN)2]C1 }'^ Although our analysis using the mean-field approximation 
of the Hubbard model may oversimplify the complexity, a more than two orders magnitude 
difference certainly supports the idea that the quantum criticality should be seen. 
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Fig. 3. Scaling of carrier density measured from critical point \X — Xc\ as functions of g — gc near 
MQCP at t' /t = 0.3645, which determines exponent 6. At r=0, MQC scaling, 6 = 2, is always 
satisfied, while a crossover from S = 2 to 5 = 3 is seen at finite temperatures. Results are obtained 
by numerically solving the Hartree-Fock equations around MQCP. 

In its literal form, the present non-GLW picture appears to apply when some order such 
as antiferromagnetism exists in the both sides of the transition. However, the opening of the 
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Fig. 4. (color): Left panel: Critical temperature plotted against gi!Tc = 0) — g at t'/i = 0.3645. 
Right panel: Ising and MQC regions separated by crossover. The crossover point is defined by the 
crossing point of two asymptotic lines in the plane oi X — Xc and \g — gc\ as in an example of 
two dashed lines in Fig. 3. The distance between the critical point and the crossover measured by 
Iff — 9c{Tc)\ (abscissa) is plotted as a function of Tg (ordinate). 



preexisting gap itself in metals is not necessarily the consequence of the symmetry broken 
order. The quasiparticle gap between the so-called lower and upper Hubbard bands may open 
without the symmetry breaking. ^^ Therefore, the present consequence may survive without 
the strict order, where A is replaced with a correlation-induced gap without long-ranged 
orders. 

When filling is controlled, the first-order MIT appears as a phase separation, ^^ and the 
MQCP appears at the boundary of this phase separation region at T = 0. Around MQCP, 
the carrier density is scaled by X oc ^^' with the chemical potential ^ and 5 = 2. This 
results in the charge compressibility k = dX/dfi oc X^^ = X^^. The compressibility also 
shows K oc \g — gc\~^ with 7 = 1 similarly to the case of bandwidth control. Therefore, k 
generically diverges in all the routes towards MQCP. It should be noted that quantum Monte 
Carlo studies^ have also shown k, oc 1/X, implying that the present result holds beyond the 
mean-field approximation. 

The present result may be justified beyond the mean-field approximation for the following 
reason: X has its intrinsic length scale ^ = 1/yX in the metallic side, because the carrier 
density X has the unique diverging length scale of the mean carrier distance, ^, yielding 
^ oc ((7 — gc)'^ oc B~'^ oc X^'^ with i' = 1/2 defined at ^ = 0. Then the dynamical exponent 
z of MQCP becomes 4 because k/X oc ^^ is satisfied.^' ^'''^^ The mean field approximation is 



8/10 



J. Phys. Soc. Jpn. Letter 

marginally justified, because the Ginzburg criterion^^ for the upper critical dimension given 
by 7 + 2/3 — i'{d + z) = is satisfied for 7 = 1, /3 = 1, i^ = 1/2, d = 2 and z = 4. The critical 
exponents then become exact except for logarithmic corrections. 

Such density fluctuations of degenerate Fermi particles may induce various instabilities 
toward orders including superconductivity, and profound effects on competing orders may 
be triggered. Our flndings open a new route of studying the physics of electron correlation 
effects. The result implies the crucial importance of enhanced dielectric fluctuations typically 
in the energy range up to the insulating gap amplitude (typically up to ~ leV) at a small 
wavenumber. Even in the fllling control transition with the constraint of the long-ranged 
Coulomb force ignored in the Hubbard model, such enhanced dielectric response may play a 
crucial role in the dramatic change of spectral properties over a large energy range as known in 
optical conductivity.^ It is desired to develop experimental probes for carrier density-density 
correlations beyond the present accuracy of available probes such as electron energy loss 
spectroscopy (EELS) and X-ray. In fact, an accurate probe should show "opalescence". 
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